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1.  Introduction 


■v 


The  distributions  of  functions  of  the  eigenvalues  of  the  real 
and  complex  Wishart  matrices  are  very  useful  in  studying  the  struc 
tures  of  the  covariance  matrices  of  the  real  and  complex  multi¬ 
variate  normal  distributions  respectively  and  other  problems. 
Krishnaiah  and  Lee  (1977)  derived  the  joint  asymptotic  distribu¬ 
tions  of  the  linear  functions  as  well  as  the  ratios  of  the  linear 
functions  of  the  roots  of  the  central  Wishart  matrix  when  the 
population  covariance  matrix  has  simple  roots.  Pujikoshi  (1978) 
derived  an  asymptotic  expression  for  the  distribution  of  a  func¬ 
tion  of  the  roots  of  the  central  Wishart  matrix  when  the  roots 
have  multiplicity  whereas  Krishnaiah  and  Lee  (1979)  obtained 
corresponding  expressions  for  the  joint  density  of  the  functions 
of  the  roots.  In  this  paper,  we  obtain  asymptotic  expressions 
for  the  joint  densities  of  various  functions  of  the  noncentral 
real  and  complex  Wishart  matrices.  These  expressions  are  in 
terms  of  multivariate  normal  density  and  multivariate  Hermite 
polynomials.  Percentage  points  of  some  test  statistics  are 
computed  by  using  the  above  asymptotic  expressions  and  these 
percentage  points  are  compared  with  the  results  obtained  by 
simulation.  Applications  of  the  above  results  are  also  discussed 
in  problems  of  studying  the  structure  of  interactions,  mixtures 
of  multivariate  normal  populations,  and  reduction  of  dimension¬ 
ality.  The  results  obtained  on  the  joint  distribution  of  the 
functions  of  the  eigenvalues  of  the  real  Wishart  matrix  are 
generalized  to  the  case  of  multivariate  quadratic  forms.  Finally, 
the  joint  asymptotic  distribution  of  the  functions  of  the  roots 

<  <  -j 

of  the  complex  Wishart  matrix  is  derived. 


1 


2.  Perturbation  Technique 


Let  ft,  _>. . . _>ft  be  the  eigenvalues  of  the  symmetric  matrix 
T:pxp,  and  A^_>. .  .>1^  are  the  eigenvalues  of  the  symmetric  matrix 
V:  pxp,  where 


T(e)  =  V  +  V 


(1) 


.2  „(2) 


(2.1) 


Then,  there  exist  orthogonal  matrices  r  and  G  such  that 
T  =  GLG '  and  V  =  fAr*,  where  L  =  diag  (ft,,..., ft  ), 

r 

A  =  diag  (A,,..., A  ).  The  columns  of  r  and  G  consist  of  the 
eigenvectors  of  V  and  T  respectively. 

Lawley  (1956),  Mallows  (1961),  Izenman  (1972)  and  Fujikoshi 
(1978)  have  approximated  the  eigenvalues  and  eigenvectors  of  T 
in  various  papers.  The  authors  have  either  assumed  that  A^fs 
do  not  have  multiplicity  or  the  approximations  were  established 
by  tacitly  assuming  that  the  eigenvalues  and  eigenvectors  admit 
series  expansions  in  the  infinitesimal  parameter  e  as  follows: 


Sj 


+  eA. 


+ 


(l) 

(l) 


+  e2A.(2) 
J 

+  eV2) 


+  •  .  . 


+  .  .  • 


(2.2) 


and  no  attempts  were  made  to  prove  that  the  series  converge. 

A  more  insight  treatment  to  settle  this  question  of  convergence 
is  found  in  Kato  (1976). 

Now  T(e)  and  V  are  linear  transformations  which  operate 
on  the  p-dimensional  complex  vector  field  Cp,  e  is  also  complex, 
Ai>. . .>Ap  are  the  eigenvalues  of  V:  pxp  such  that 
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ior  a  =  ±, r  ,  q^+...+qr  =  p  ,  qQ  =  0  and  let 
J  (a=l,...,r)  denote  the  set  of  integers 

q^+. . .+qa_1+l ,  . ..,  q1+...+qa.  We  need  the  following  lemma 

in  the  sequel. 

Lemma  2.1.  For  Hermitian  matrices  T(e)  and  V  as  defined 


before. 


T(e)  =  V  +  e  V(1)  +  e2  V(2)+. . . 


and  V  is  diagonalized  as  V  =  diag  (X1,...,Xp).  Then  the 
mean  eigenvalue  of  T(e)  corresponding  to  0  which  is  the 

(X 

eigenvalue  of  V  with  multiplicity  qa  ,  is 

*a*(<=)  =  0a  +  e  5Si1)  +  e2  l[2)  +‘*-  ( 2.* 

where 


j(D 

a 


j(2) 

a 


~  tr  V(1) 
Qa  aa 


=  ~  tr[V^2)  +  l  (I 

n  rv  rv  v 


a  x-1  v(l)  v(l)n 
%  >  Va3  V3  a  3 


(2.5: 


0-0. 
a  3 


with 


V  (  i ) 

11  v12 


r(i)  _ 


v(i>  VU) 

„  rl  vr2 


.  V  ^  I 

rr  J 


3 


and  ls  of  order  QaxQg* 

When  q1=...qr=l,  the  above  lemma  was  proved  in  Lawley  (1956). 
When  qa>l  (ot=l,...,r)  ,  the  lemma  was  given  implicitly  in  Kato 
(1976).  For  qQ=l  the  normalized  eigenvector  of  TGe),  .corresponding 
to  6a  is  Ga(e)  =  (Gla(e) , . . .  ,Gpa(e.) )  ’ ,  with 


a,  (e) 

e(e)  =Te  "-A  )  +  I 

a  i  j^a 


aij(E)aJa(e) 


( e  -  a  . : 

a  1 


e  -a, 

a  j 


alc,U)aac,(e) 

<W2 


. .  .i^a 


1  ”  2 


j  aim(e)  aml(e> 

( W(  V*? 


(2.6) 


where 


A(e)  =  T( e)  -  V  =  (a^U)) 

The  series  (2.4),  (2.6)  are  convergent  for 

2c,  , 

|e|  <  (—  +  c 2)~x  (2.7) 

where  c1,c2  >  0  such  that  ||V^||  <  c-jC^”1  for 
j  =  1,2,...  ,  and  d  =  min(|6a  -  ,  1 6a  -  Qa+1\). 


•i ff 

3.  Asymptotic  Joint  Distributions  of  Functions  of  the 
Roots  of  Noncentral  Wishart  Matrix 

Let  the  columns  of  X:  pxn  be  distributed  as  multivariate 

normal  with  covariance  matrix  E  =  (c.  .)  and  means  given  by 

t  J 

E(X)  =  U  =  where  uj  =  (p^ ,  •  •  • ,  pJp) .  Then, 

S  =  XX'  =  (S.  .)  is  distributed  as  the  central  or  noncentral 
ij 

Wishart  matrix  W  (n,E,M)  with  n  degrees  of  freedom  according 

P  n 

as  M  =  0  or  M  f  0  where  M  =  l  p  .p.'.  =  n(v.  . ) .  Now,  let 

~3~J 

_>•••_>  denote  the  eigenvalues  of  S/n  whereas 
denote  the  eigenvalues  of  E(S/n)  =  E  +  M/n  =  A  .  Without 
loss  of  generality,  we  assume  that  A  =  diag.  ( A ^ , . . . , Ap) . 
Also,  let 


V 


=:•••=  X  =  B 

■+q  ,+1  q.+**»+q  a 

Ma-1  M1  a 


(3.1) 


for  a  =  l,2,»**,r,  q1  +  •••  +  qr  =  p,  and  qQ  =  0. 

In  this  section,  we  derive  the  joint  asymptotic  dis¬ 
tribution  of  ,  •  •  • ,  L^  where  L..  =  /n  {T^  (  5^  ,  •  •  • ,  2.p)-Tj(  A^  ,  . .  .  ,  Ap) } 

and  T £  )  satisfy  the  following  assumptions: 

J  A  P 


(i) 

Tj(V  ***’ V  iS 

analytic  about  A^, 

(ii) 

dTi(A1> •••,ap) 

dl. 

J1 

S2Tt( 

k=A  ijl  ia 

1 

3?,.  3*.  I  =  C  i  3 1 J  2  aia3 

H  31  l=\ 


5 


cr 


3' W“**V 

3 1.  32.  .  32.  - 


2=  A 


Cij1J2j3  aia3y 


(3.2) 


for  j1  e  Ja,  j2  e  Jg,  J3  e  J^A'  =  ( Xi » *  *  * » Ap> * £ '  *  (4i»***’4p> 
and  denotes  the  set  of  integers  q1  +***+Qa_^  +  1 > *  *  * >^1 +  *  * *+  Q 


a 


a 


for  a  =  1,2, • » • ,r. 

Expanding  Tj[(  >  “  * »  *p)  as  the  Taylor  series,  we  obtain 

VV-’V  -'TiOr—-V  +  <V  V 

1  Gt 


.  r  r 

+  \  l  l  aiae 

z  ct=l  8=1  lfltp 


I 


l  <*s  -68> 


3leJB  VJB  'h  “  h 


,L  3L(V“,(W(VV 


aBy 


j Jr>  J2tu8  j 3 t-w y 


+  terms  of  higher  degree. 


Now,  let 


where 


S/n  =  A  +  -  V 

/n 


(3.3) 


(3.4) 


V  =  /n  (|  -  A)  * 


,V11  V12 


lr 


V21  V22  V2r 


\ 


\  *  *  ’  / 
'  V  .  V  ~  •  •  •  v  / 


Vrl  Vr2 


rr 


So  by  applying  Lemma  2.1  on  (3.4),  we  obtain 
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XU  ^ 


✓n  a=l  iU' 


r  r 

l  l 

o=l  3=1 


(tr  Z^Ktr  Z(1))>  +  0(n-1)  (3.5 


where  Z(1)  =  V  ,  Z^2)  =  T  0“i  V  0  V„  and  0  , 
a  aa’  a  a3  a3  0a  a3 

A]  so. 


0  -  0Q, 

a  3 


tr  *<*>  -  rff  l  (-2i=i  -X.  V 

j,  cJ  \  n  J1  / 

1  a 


tr  z(2)  =  -  J  e“*  l  Is2 

a  n  a3  •  f:T  •  j;  t  J- 


3*a  jl£Ja  j2eJ3  J1J2 


After  some  algebraic  manipulations,  we  obtain  the  following 
expression  for  the  joint  characteristic  function  of  ,  •  • « , : 


'(tl'***,tk)  =  E{exP(i  3^  tjLj)} 

=  Efexpd  I  l  t  a  tr  z'15) 

j=l  a=l  J  J 


(3.6 


k  r 


k  r  r 


U+  -Hi  l  t  a  tr  Z<2)  +|  l  l  It 
✓n  j=l  o=l  J  Ja  a  2  j=l  o=l  3=1 


tr  tr  z£X)) 

a  p 


+  0(n_1) } 3  =  E1(t)  +  E2(t)  +  E3(t)  +  OCn”1) 


where  t'  =  (t^,»»«,t^).  In  Eq.  (3.6), 


k  r 


(Dm 


!,( t )  =  ECexp( i  l  l  tiaia  tr  za  )] 
1  ~  i=l  a=l 


(3.7) 


where  B 

h<V 


E3(t)  - 


=  etr  (-i  /n  B  )  [  I-2i  B  £/ /n  |  ^ 


x  exp(i  tr CM( I-2i  B  £/ /n)  ^  B  //n  )} 


l  ti  diag(ci;L>  •  •  •  ,cip)  •  Also, 


k  r 
E[ 

/n 


k  r 


Z  l  l,  tiaia  tr  Zi2)  *  eXpU  ihah  Wal 
n  i=l  a=l  11 


tr 


a. 


=  Ej.  (t)  -  l  l  l  l  l  tjaia  °aB 
/n  i=l  ot-1 


x  1  r  v  /a*  a*  +a*2 .  +a*  -  £?■  •  +2aT  . 

n  -jlx  3^1  J2j2  J1J2  J1J1  332  J1J2  jjl  jj2 


*  r-2  . 

+  a  .  •  5  .  .  ) 
J2J2  JJ1 


+  cjIi<aJ^2+eJJieJ^>l2> 

EC^  !  I  1  VlaB<tr  Zi1>)(tr  *»l>> 

2/n  i=l  a=l  3=1 


(3.8) 


k  r 


expU  l  l  ti  ai  a, 
i^=l  a1=l  111 


tr 


a. 


E  (t)  —  III.  I  .1,  tiaia8 

1  ~  2/n  i=l  a=l  3=1  3ieJa‘,#J 8 


{I  ?  (2a*2i  +  4ot  .  £ j j  ) 

^  3 1^2  J1J2  Jjl  332 


8 


) 


n  n 


4.  i.  V  V  fn*  4.r*  *  9  n  *  2 

n  j=l  m=l  h*k  jjl)(oj2j2Hmj2)  "  XjuI1(0j2j2Hjj2) 


-  X.  I  (o*  .  H2..  )  +  n  X.  X.  } 
J2  j=j_  JJj  J1  J2 


(3.9) 


where 


* 

E  =  E 


2iB1E 

-1 

\  _ 

(o*  .) 
iJ 

/n 

) 

2iEB1  ^ 
/n  J 

-1 

1 

ti 

C-i. 

»-* 

(3.10) 


Using  the  expansion  that 


|  I-A |  ^  =  exp  0(  l 


tr  A^ 


3=1  3 


in  (3.7)  and 


(I-A)"1  =  l  AJ 
j=0 


in  (3.10),  we  obtain 

k  p 

o*  =  o.  .  +  —  Y  f  t.c..  a.  .a.  .  +  0(n_1) 

jlJ2  J 1^2  /n  i=l  j=l  1  1J  J1J  J2J 

F  .  =  u  .  +  —  y  f  t.c.  a.  p.  +  0(n_1) 

JJl  «1  /n  i=l  m=l  1  lm  Jlra 


Eq.  (3.7),  (3.8),  (3.9)  lead  to 


f(t)  =  exp( -  |  t ’ Qt ) 


u  *ij  i  ti(d1+d2)  ♦  i  !  ti3ti  * ,  « i/ws*1 

/n  i=l  1i^2i3  1  Z  ° 


+  0(n-1 ) 


(3.11) 


fil>  (i2)  (i,)  (i9) 

where  Q  =  (Q.  i  ) ,  Q.  •  =  2tr  R  R  +  4tr  R  1  ill  z 

ll  2  X1  2  i 

and  Q  is  assumed  to  be  nonsingular .  Also, 


r  r 


dl  “  l  1  l  I  aict  °afi(oi  i  ai  i  +ai  i  +2oi  i  v-i  -i 

a=l  Bfa  jjeJa  j2cJ$  K  J1J1  J2^2  J1J2  J1J2  ^1^2 


a .  .  v .  .  +a .  .  v .  .  ) 
J1J1  J2J2  J2J2 


r  r 


d  =  y  y  y  y  a .  „(o.  .  +20.  .  v.  .  ) 


«=1  8=1  jl£Ja  j2cJ6 


Jl‘]2  ^1^2  J1J2 


4  „(l2}  J13) 


(il)  „(i2>  /V 


gj  =  |  tr  R  x  R  R  +  4  tr  R  x  R  if,  ° 

j.  ( i2)  ^  -i2)  ^  ^2^ 

e2  “  4  I  I  I  I  ai  a  9a8(5i  i  +Ti  i  +Ti  i  ^ 

a=l  B^a  j  eJ  j0cJQ  1l“  afS  J1J2  J1J2  J2J1 

1  Ct  &  p 


(io)  (i„)  ( io) 

x  (5  d  +T  J  +T  *  ) 

M2  M2  Mi 


r  r  ( i9)  (i9)  (io)  (io) 

2  I  I  .  I  I  a.  (E_  +2T  •  ) ( H  +2T_  ) 

g  1  6-1  j1eJa  j2eJ6  1  ‘Ml  Ml  ^2^2  M2 


(3.12) 


We  define  here  M/n  -  (v  ) 

=  diag(cil , • • • ,cip) 

r(  i)  =  E,  n, 


EC(iM  T(l) 


=  -M  1  ^ 

n 

(3.13) 


where  A.  .  denotes  the  (i,j)th  element  of  matrix  A  =  (A..) 
1  j  ■*■  J 
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*<A  <  "  ‘ 


Now  inverting  (3.11)  we  obtain  the  following  expansion 
for  density  of  L  =  (L.  ,  •  •  • ,  L  ) 

1  K 

f < L1 , • • • , Lk )  =  N(L,Q)  x 
1  k 

[1  +  ~  l  i  H  (L)(d  +d  ) 

/n  i=l  1  1  2 

+  ~Z  f  Hi  i  j  (LMg^+go+g^.))  +  0(n-1)  (3.14) 
/n  i1.i2*i3  1  27'3  x  2  6  ’ 


wv  3re  N(L,Q)  =  - jrri - r™  exp(-  £  L'Q_1L) 

(2n)k/2|Q|1'2  2  ~ 


H. 


.  ...  .  (L)N(L.Q)  =  (-1) 

Jl’  ,Js  ~ 


s 


(3.15) 


9L .  9L. 

V*’  Js 


N(L,Q) 


Now,  let  TiU1»,##»*p)  =  A±.  Then  L±  =  /n  (J^-Xj.). 

Using  Eq.  (3.14),  we  obtain  the  following  expression  for  the 
joint  density  of  the  roots  J^,***,*  when  l  =  o2I : 

f(Li»,,**Ln)  =  n^l.QH1  +~  f  H.(L)  l  e:J(o2x,+o2x,-o4) 

p  ~  /n  i=l  1  ~  j^i  1  -1 

+  —  !  H  (L)(X.-  §  o2)}  +  0(n-1) 

✓n  i=l  111  "  13 

(3.16) 

where  Q  =  diag(Q^ , • • • ,Q  )  and  =  2o2(2Xi~o^).  When  £»cr2I  and 

2 

X1>. . . >Xt =  Xt+i=  ...X  -  o  ,  Carter  and  Srivastava  (1979) 

obtained  an  alternative  expression  for  the  joint  density  of 
£p  by  using  a  different  method. 


11 


t 


The  general  k-dimensional  Hermite  polynomial  of  order  s^O 

is  denoted  by  H.  ,  i  d;A)  where  ~  = 

11» 12  ’  ' *  * »  s 

is  the  polynomial  variate  and  A  =  (fi|j)  a  positive- 
definite  kxk  matrix,  0  <  i^  <_  k  for  j=l,...,s  and 


H 


» *  •  •  * 


(x,A)  N(x,Q)  =  (-1)° 


‘  8 


3x.  . . . 3x . 
X1  xs 


N(x ,Q)  (3.17) 


where 


N(-’Q)  •  (21,  6XP(~*  -’9-1  ?> 


(.3,18) 


and 


A  =  Q 


-1 


2  J. 

For  dimension  k=l ,  Q=x"  is  a  scalar  and  A=6=  -rr 

T 


H1(x,6)  =  xA 

Hli;1(x,A)  =  x363  -  3  x  52 


(3.19) 


and 


xa 


H1(x,A)  N(x ,Q)dx 


=  -  i  x  exp  Fir) 


(3.20) 


a 

/  —  CO 


Hm(x,A)  N(x  ,Q)dx  =  d-a2)  exP  FH 


rl  Pt1t2^ 

For  the  dimension  k=2,  let  Q  =  2 

pTlT2  t2 


Then 


A  "  Q"1  “  TTT  72 


\-pTlT2 


pTlT2) 

2 

|  S3 

[«11 

6X2 1 

X1  i 

i  621 

622i 

K3.21) 


12 


Also,  we  have 


,A)=Xi  6U 

+  x2612 

(x, A)  =  x3 

6I1+  3  : 

3 

.3 

+  x2 

612~  3  ' 

3 

2 

)  r-i 
X 

II 

<3 

X? 

511  512 

6‘J,+  3  x^x25ix<S12  +  3  xi  xl2fill5!9 


.2 

:1  611' 


k2  °11  12 

.  .2  .,2 


+  x1x2(612  +  25xx5i2l522^  +  x2612  ^22 


"  3  X1  611  512  “  x2(2512+iS11  622)‘ 


(3.22) 


Similar  equations  for  H2(x,A)  ,  H222(x,A)  and  Hx22^x,A^ 
are  obtained  by  interchanging  subscript  1 ,  2  in  H^(x,A), 
H1;11(x,A)  and  Hn2(x,A)  respectively. 


Nov/ 


Tob 


T,a 


1  x*'  x®  N(x ,Q)dx^dx2  = 


T1  X1  T2  x2  N(x>^^dxidx2 


—  CO  —00 


by  changing  of  variables  and  R  = 

matrix. 

Define 


1  P 
P  1 


(3.23) 


is  the  correlation 


U 


r ,  s 


a 


x^  x2  N(x ,  R)dx1dx2 


(3.24) 


Integrating  by  parts,  we  have 

y.  n  =  -  ($(a)*(B)  +  p«Kb)*(A)) 

*•  f  ”  (  o  « ) 

p  =  — pt>4>(b)  <J>  (A)  -  pa<t>(a)  <I>(B)  +  (1-p2)  N(a,b;R) 

X  f  1 
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BJE 


]i2,0  =  U0  q  -  a  $  (a)  4>(B)  -  p  b  4>( b)  <KA) 

+  p(l-p2)  N(a,b;R) 

In  general,  the  recursive  relation  is 

"r,s+l  “  b  ’’p.s  ‘  Kr.s-1  +  •V+i.s  p 

(1-p2)  b(s-l.)ur  s  2  -  bppr+l  s  '  ( 

pr+l, s  -  a  »r,S  =  (1-p2)l,Vl,S+  Pr,s+1  p  -d-p2)a  (?-D 
*  11  i-2,s  -a  p  pr-l,s  +x 


-Ca2+2)  0(a)  0(B)  +  p( p2-p2b2-3)  0(b)  0(A) 
+  (a+pb)  p(l-p2)  N(a,b;R) 


u2,l  -  "p(a  +2>  ♦(»)  0(B)  -  (1  +  p  +  p2  b2)  0(b)  0(A) 

+  (l-p2)(a+pb)  N(a,b;R) 


vhere 


A  =  (l-p*)-2  (a-pb)  ,  B  =  (l-p2)"^  (b-pa) 


4>(a)  =  — exp  — 

/21t  2 


♦  (A) 


a 

» 

«  — C 


4>(t)dt 


^o,i*  y0,2’  ^0,3  an^  2  are  obtained  by  interchanging  a  with 
b  and  A  with  B  in  n,  y0  n,  y_  n  and  y0  .  respectively, 


4.  Applications  in  Investigation  of  the  Structures  of 

Interactions 


In  this  section,  we  discuss  some  applications  of  the 


results  of  Section  3  in  studying  the  power  functions  of 


various  tests  for  the  hypothesis  of  no  interaction  in  two- 


way  classification  model  with  one  observation  per  cell. 


Consider  the  model 


x..=  u+a.  +$.  +n-.+e-. 
ij  13  ij  ij 


(4.1) 


for  i  =  l,...,u,  j  =  1,2,. . . , s .  Here  x . .  denotes  the 

1 3 


observation  in  i-th  row  and  j-th  column,  u  is  the  general 


mean,  a.  denotes  the  effect  due  to  i-th  row,  ft.  denotes  the 
^  3 


effect  due  to  j-th  column  and  n.  •  denotes  the  interaction  of 

13 


i-th  row  and  j-th  column.  Also,  we  assume  that  e. ,'s  are 

t  J 


distributed  independently  and  normally  with  mean  0  and 


variance  o 


s 

=  l 

j  = 


Now,  let  d .  .  =  x .  .  -  x.  -  x.  .  +  x..  .where  sx.  .= 

««  ***  J  J  *^11  O  *  **• 


l  X.,,  UX  =  Six.,  and  usx..  =  l  I  x 
1  =  1  XJ  3  i=l1J  i=l  j=l  13 


Also,  let 


D  =  (d. .),  X  =  (x. . ) ,  W  =  C'  X  C  C '  X'  C  where  C  is  chosen 

1 J  1 J  USS  U  U 

such  that  Ci'jCu  »  Iu_1  and  CuC^  =  Iu  -  ^  Ju  where  Ju  is  the 

u xu  matrix  with  all  its  elements  equal  to  unity.  The  non- 


u  ^u 


zero  eigenvalues  of  DD  are  the  same  (e.g.,  see  Johnson,  and 


Graybill  (1972))  as  the  nonzero  eigenvalues  of  W.  Also,  the 


columns  of  CuX  are  distributed  independently  as  multivariate 


normal  with  mean  vectors  given  by  E(CuX)  =  CuMQ  and  a 


common  covariance  matrix  0  Iu_j  where  MQ  =  (m^j)  anc*  mfj 


=  n  +  6j  +  nj.1 


In  addition, 


mam 


7 


E(W/(s-l)) 


2 

o  I  ,  + 
u— 1 


I  t 


_  r  (4.2) 

**0 

and  C  M  C  C  M  C  =  C  nn  C  =  ft  where  n  =  (n. •).  So,  W  is 
u  o  s  s  o  u  u  u  ij  ’ 

distributed  as  the  noncentral  Wishart  matrix  with  (s-1) 
degrees  of  freedom  and  noncentrality  matrix  ft.  When  n=0, 

W  is  distributed  as  the  central  Wishart  matrix  with  (s-1) 
degrees  of  freedom. 

Let  1  >  •  •  •_>  £u_1  be  the  eigenvalues  of  W/(s~l)  and  let 
be  the  nonzero  roots  of  £  .  Then,  the  problem 
of  testing  the  hypothesis  H:  ft  =  0  is  equivalent  to 
testing  the  hypothesis  that  the  eigenvalues  of  are  equal. 
Suppose  n  =  X  ct0 '  where  a1  =  (a.,..., a  )  and  3'  =  (3,,..-,B  )• 

~  ~  ~  -L  l\  ~  1  S 

2  1 

Then  ft  =  A  C  ag'Ba'C  ,  and  the  nonzero  root  of  ft  is 
u~~  ~~  u 

A3 '3  a 'a.  Next,  we  will  assume  that  the  rank  of  n  is  c. 

Then,  using  the  well-known  singular  value  decomposition  of 
the  matrix,  we  can  write  n  as 


A  .  w1v1  +  .  .  •  + A  w  v 

j.  ~  .1  ~  1  (!  ~  C  ~  C 


(4,3) 


2  2 

The  nonzero  eigenvalues  of  nn '  are  A ^ ,  .  .  .  ,A  and  the  asso¬ 
ciated  eigenvectors  are  w.  ,  .  .  .  ,  w  •  The  eigenvectors  of 

~  X  ~  o 

2  2 

n'n  corresponding  to  the  eigenvalues  A ^  ,..., Ac  are  v^,...,^. 

2  2 

The  nonzero  eigenvalues  of  ft  are  \1,...,A(,. 

The  problem  of  testing  the  hypothesis  of  no  inter- 

I 

action  in  two-way  classification  with  one  observation  per 
cell  was  studied  by  several  authors  (e.g.,  see  Tukey  (1949) 


1 6 


and  Williams  (1951)).  The  statistic  proposed  by  Tukey 


^  A  A 


is  given  by  (a'nB)  /(a'aKB'B)  where  a’  =  ), 

$'  =  (S1,...,Bs),  n  =  <  n±j ) ,  ou  =  x.>  -  x..,  §  =  x..  -  x.. 

A  _  _ 

and  =  x.^  -  x^  -  x . ^  +  x..  .  Gollob  (1968)  and  Mandel 
(1969)  considered  the  problem  of  testing  the  hypotheses 
A.  =  0  individually  under  the  model  (4  .3)  by  using  the 

J 

statistics  F.  =  Z./  (£.+...+£  ) .  Gollob  treated 

J  J  U“*l 

i  as  independent  chi-square  variables  to  get  an 
approximation  to  the  distribution  of  F  ^  .  But  ’ s  are 
neither  independent  nor  distributed  as  chi-square  variables. 
Corsten  and  Van  Eijnsbergen  (1972)  showed  that  the  like¬ 
lihood  ratio  statistic  for  testing  the  hypothesis 
A^=...  =  Ac  =  0  under  the  model  (4.3)  is  (£^+...+J,  )/ 

(Z  .).  When  c=.l ,  this  statistic  was  derived 

c+1  u-1 

independently  by  Johnson  and  Graybill  (1972).  Schuurmann, 
Krishnaiah  and  Chattopadhyay  (1973)  and  Krishnaiah  and 
Schuurmann  (1974)  discussed  the  problem  of  testing  the 
hypotheses  A^  =  0  simultaneously  by  applying  the  simultan¬ 
eous  tests  of  Krishnaiah  and  Waikar  (1971)  for  the  equality 
of  the  eigenvalues  of  the  covariance  matrix  of  the  multi¬ 
variate  normal  population.  Ghosh  and  Sharma  (1963)  studied 
the  power  function  of  Tukey 's  test  for  n^.  =  0  against  the 
alternative  that  =  A  ft  ^ .  Yochmowitz  and  Cornell 
(1978)  derived  the  likelihood  ratio  test  for  A^  =  •••  =Aa  =  0 
against  the  alternative  that  Aa  f  0  and  Aa+-j=  •  . -=AC  =  0. 

We  now  compare  the  power  functions  of  various  procedures  for 
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T 


testing  the  hypothesis  of  no  interaction. 

Let  Ti=  t1/au_J,  T2=  (tr  W/u-J )U_J / | W | ,  T3=  (u-c-1)  * 

£l/(£c+l+' ' ,+£u-l) ’  T4=  lj.+  -  •  -  +  ^c)/c(  JLc+i+*  •  *+ilu_i)  • 

2 

When  o'  is  unknown,  we  can  use  any  of  the  above  statistics  for 
testing  the  hypothesis  of  no  interaction. 

If  we  use  T  we  accept  of  reject  according  as 


T.  >  c  (4.4) 

i  a 

where 

PCTi  -  calH0]  =  (JL"a)  (4-5) 

The  test  statistics  is  based  upon  the  statistic  con¬ 
sidered  by  Krishnaiah  and  Waikar  (1971)  for  testing  the 
sphericity,  whereas  the  test  statistic  is  bused  upon 
the  likelihood  ratio  test  statistic  for  sphericity.  The 
statistic  j.s  the  likelihood  ratio  test  statistic  (see 
Corsten  and  Van  Eijnsbergen  (1972))  for  testin'  „e  hypothesis 
of  no  interaction  of  multiplicative  compone’  model  (4.3). 

Table  1  gives  a  comparison  of  the  power  functions  of 

various  procedures  for  testing  the  hypothesis  of  no  interaction 
2 

when  o  is  known.  The  rows  corresponding  to  S  denote  the 
simulated  values.  The  multivariate  normal  deviates  are 
generated  by  the  IMSL  subroutine  GGNRM,  and  10,000  trials  are 
performed  for  each  case,  the  95%  confidence  limit  for  each 

/V  X  A 

value  is  then  1.96{£(l-p)/l0,000}“ ,  where  p  is  the  actural 
value  from  the  empiricaJ  trials. 


18 


in 

The  rows  corresponding  to  N  denote  the  values  corresponding  to  the 
first  term  in  the  asymptotic  expansion.  The  rows  corresponding 
to  N  +  0(n_i)  give  the  values  corresponding  to  the  fust  two 
terms  of  the  expansion. 
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The  radius  of  convergence  is 


where 


?\  |v| 


V  =  /n  (S/n-M) 


(4.6) 


We  will  choose  n  such  that  l//iT  K  Yq  Now  j  ( V  j  j  =  CL( /n( S/n-M) ) , 

where  CL(  • )  denotes  the  largest  root  of  /n  ( (S/n)-Nl)>  and  is 

2  2 

approximately  distributed  with  mean  O.and  variance  2o  (X^-o  ). 

In  Table  1,  consider  the  entry  when  p=3,  XA=12, 

2 

X_  =  6  and  X0  =  1.  In  this  case,  o  =  1,  d  =  XP  -  Xq  =  ^ , 


with  one  standard  deviation, 

d 


2v2cr2(2A1-a2) 


/7TS 


where  n=10  is  chosen. 


When  the  entry  is  p  =  4  ,  X^  »  2.5f  \^~  1.7  and  Xg  =  1, 

then  02=  1,  Xj,  =  2.5,  d  =  X^  -  X^  =  0.8  and 
n  =  50. 

The  table  reveals  that  results  based  on  normal  approx¬ 
imations  are  not  sufficiently  accurate  for  n  as  large  as 
100,  while  the  asymptotic  expression  taking  the  term  of 
order  n~~  achieves  numerical  accuracy  for  moderate  sample 
sizes.  This  suggests  that  care  should  be  given  for  the 

statistical  inferences  which  are  based  on  the  normal  approximations. 

2 

Next,  consider  the  model  (4.1)  when  a  is  known.  In  this 
case,  we  accept  or  reject  the  hypothesis  X^=...=Xc  =  0 


according  as 


(4.7) 


where 


p  t  ~j  1  d2a|H]  =  (1-a)  (4. 

o 

When  H  is  rejected,  the  hypothesis  0.  =  0  is  accepted  or 

2  < 

rejected  according  as  (JL/o  )  >  d2o[.  When  H  is  true,  is 
the  largest  eigenvalue  of  the  central  Wishart  matrix.  Exact 
distribution  of  this  statistic  is  given  in  Krishnaiah  and 
Chang  (1971)  and  exact  percentage  points  are  given  in 
Krishnaiah  (1980).  When  II  is  not  true,  an  asymptotic  expres¬ 
sion  for  the  distribution  of  l ^  can  be  obtained  as  a  special 

case  of  (3.14)  if  A,  is  different  from  0~,...,0  . 

i  "  r 

When  A1>...>A2>0,  Srivastava  and  Carter  (1980):  obtained 
asymptotic  expression  of  log( £./£.+... +£  )  andT^^"1^ 

by  a  different  method.  For  a  review  of  the  literature  on 
tests  for  no  interaction  in  two  way  classification  model 
with  one  observation  per  cell,  the  reader  is  referred  to 
Krishnaiah  and  Yochmowitz  (1980). 


5.  Applications  in  Cluster  Analysis  and  Reduction 

of  Dimensionality 

LcL  be  Independent.  p-d  itnensiotr.il  random 

variables.  We  consider  using  the  sample  covariance  matrix 

N 

s  =  l  (x.-X)(x.-x) • 

i=l  ~  -1  - 

where  X  »  N~*(X^+. . ,+X^) .  We  wish  to  test  the  hypothesis 
that  X^'s  come  from  a  single  multivariate  normal  population 
with  covariance  I  against  the  hypothesis  that  they  come 

from  a  mixture  of  k  £  p  such  populations,  differing  in 

2* 

means.  We  assume  T.  =  a  I.  For  k=2  the  null  hypothesis 
and  the  alternative  hypothesis  Hg  are  given  by 

Hj  :  X.  -N(p,02I) 

H2:  X±  -  wN( iij , a2I )  +  (1-ir)  N(]jg,a2I) 

where  it  is  the  mixing  probability.  Under  II?  it  is  known 

(e.g.  see  Bryant  (1975))  that 

N 

S~  l  (h  TrJ(l-77)N_jW,  (N-l;cr2I,M,) 
j=0  J  p  J 

Mj  ==  N-^ j (N- j  )  (  Uj-Rg)  (  Hl“p2^  ' 

2  2 

M.  is  of  rank  1.  Now,  let  A  =  (y  -u,,) '  (y-, -Vo)/o  which 
J  ~1 

is  proportional  to  the  largest  root  of  Mj/(N-1).  When  the 
null  hypothesis  is  true,  we  know  that 

S  ~  W  ( N-l , 02J , 0) . 

Let  the  test  statistics  T1  and  r2  be  given  as  below: 


*  "V  '^'"^^y'a  «&n  ?&*?$$&  *^9^385^5^ 


Tx  =  ^/a 


(5.1) 


T2  = 


(p-DJtj 

S,2+“  ,  +  ?p 


where  S,^_>. .  .>£  are  eigenvalues  of  S/(N-1). 

If  we  use  the  statistics  T.  Lo  test  ,  then  we  accept  or 


reject  according  as 


T-  <  c  . 
i  >  oti 


(5.2) 


(5.3) 


where 


P{  Ti  i  cai  l”l}  ‘ 


(5.4) 


Let  f^(*)  be  the  asymptotic  density  of  a  function  of  the 
eigenvalues  of  S  /(N-J.)  when  j  of  the  samples  come  from  population 

1.  Under  this  condition  ,$  ~  W  (N-1,cj2I  M.).  So 

i  P  ’  ’  j ' 

the  unconditional  asymptotic  density  function  of  the  function 
of  eigenvalue  of  S  is 

N 


?  ("I  (.(•) 

1=0  11  J 


The  following  table  gives  a  comparison  of  the 
asymptotic  power  value  with  the  simulated  value  of  tests 
of  against  Hg  for  r=4 ,  <x=0.05,  tt  =  mixing  probability, 
A  »  IIp-l-I^H/0  •  N  =  samPle  size. 

it  =  .25  it  *=  .50 


Test 


A=1 


A=1 


Tl 


M=  51 


Simu . 


.02  .42  .95 
.07  .45  .94 


.02  .63  1.00 

.07  .62  1.00 
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it  =  .25  it  =  .  50 


Test 

A=1 

2 

3 

A=1 

2 

3 

t2 

.07 

.68 

.99 

.11 

.88 

1.00 

Simu. 

.11 

.70 

.98 

.14 

.86 

1.00 

Tl 

.05 

.76 

1.00 

.08 

.96 

1.00 

Simu. 

.13 

.78 

1.00 

.21 

.96 

1.00 

T2 

.13 

.93 

1,00 

.22 

1.00 

1.00 

Simu. 

.13 

,93 

o 

o 

.27 

1.00 

1.00 

When  A  =  1,  the  largest  roots  of  Mj/(N-1)  are  clo.se  to 
zero,  the  radius  of  convergence  for  perturbation  approx¬ 
imation  of  eigenvalues  is  small,  and  the  asypmtotic  expression 
does  not  give  good  approximation.  Note,  that  if  the  means 
under  H2  are  separated  by  more  than  two  or'  three  standard 
derivations, that  is, for  A  =  2,3,  one  may  reasonably  expect  to 
detect  the  presence  of  two  components,  while  if  they  are 
separated  by  less  than  two  standard  derivations  the  detection 
generally  will  not  be  gooh. 

Consider  k  p-variate  normal  populations  with  unknown 
mean  vectors  VU  , . . . ,  p,,  and  a  common  known  covariance  matrix 
E.  We  assume  that  n^  observations  are  available  from  i-th 
population  and  the  sample  mean  vector  and  corrected  sum  of 
squares  and  cross-products  (SP)  matrix  based  on  these  obser¬ 
vations  are  respectively  given  by  x^  and  .  The  SP  matrix 
which  explain  the  variation  between  groups  is  given  by 
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"l?i. 


I  •  =  I  I  x. 

1=1  3:1  1=1  .1=1  ' 


(5.5) 


and  n  =  n^+.-.+n^.  Now,  let  •  •  >£  den°t°  the  eigen¬ 

values  of  BE-*.  Fisher  (1939)  studied  the  problem  of  whether 
there  are  any  structural  relations  between  the  p  components 
of  mean  vectors.  This  is  equivalent  to  testing  the  hypo¬ 


thesis  Hq  where 


:  Cy, 


(5.6) 


for  i  =  1, k  where  C:  sxp  and  £  are  unknown  and  the 
rank  of  C  is  s.  We  assume  that  t <k-l  where  t  =  p-s  The 
likelihood  ratio  statistic  for  testing  Hq  is  given  by 
Ui  =  (^t+i4, • ,+^p^ *  A  detailed  discussion  of  the  above  pro¬ 
cedure  was  given  in  Rao  (1965).  The  statistic  BE-1 
is  distributed  as  the  noncentral  Wishart  matrix  with  y  degrees 


of  freedom  and 


where  Y  **  k-1 


E(BE  *)  =  (k-1)  I  +  flE-1  =  E. 


=  I  n. (y.-y)(y  .-y. ) ’ 

X- *1  x  ~  j.  •** 


(5.7) 


(5.8) 


and  ny  =  n^y^+.  .  ,+n^y^.  Also,  let  X1  >_•••>!  denote  the  roots 
of  E*.  The  distribution  of  can  be  obtained  as  a  special 

case  of  (3.14)  and  so  the  power  of  the  likelihood  ratio  test 
for  Hq  can  be  studied  using  the  results  in  this  p^per. 


6.  Asymptotic  Joint  Distribution  of  Functions  of  the 
Eignevaiues  of  Multivariate  Quadratic  Form 


In  this  section  wo  shall  derive  the  joint  asymptotic 
distributions  of  functions  of  eigenvalues  of  the  multivariate 
quadratic  form  S  =  XGX' ,  where  we  assume  that  Gs=  0(1),  for 
whatever  power  s  raised  on  G.  G  is  a  symmetric  matrix  and 
the  columns  of  X:  pxn  are  distributed  independently  as 
multivariate  normal  with  covariance  matrix  E  =  (cKj)  and  means 
E(X)  =  U  =  (v^ . yn)  .  Then 


E(s/n)  =  ~  T.  +  =  A  =  diag(Xr...,Xp), 


(.6,1) 


Also  assume  that 


\  +  -..+ Vi+1  +  e“ 

for  a=l,2,...,r,  q^+...+q^=  p,  and  q0=0  and  let  Jl^.  •  ._>£ 
be  the  eigenvalues  of  S/n. 

We  consider  the  joint  asymptotic  distribution  of 

where  Lj  =  /n  (Tj  (  ,  .  .  .  ,  £p)  -  Tj  (  X1 . Xp) ) , 

which  satisfy  assumptions  of  (3.2).  Let 


(6.2) 


S/n  =  A  +  V//n 


(6.3) 


Using  the  Taylor  expansion  of  T^(  ,  .  .  .  ,  £  )  for  i=l,2,...,k, 

as  in  (3.3)  and  the  application  of  Lemma  2.1  we  obtain  the 
same  equations  as  (3.5),  (3.6)  and  the  characteristic  function 


nt.,...,tk)  =  E1(t)  +  E2(t)  +  E3(t)  +  0(n-i) 
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E1(t)  =  exp(-i  /n  tr  I^A)  |l  -  2i(G  (g)  B^/Zn)  |“* 

expTi  tr  C5'(I  -  2i  G  (g)  B^/Zn)-1  G®B1//nl 


i  k  r 


V9  =  E1{V  SZ  l  l  l  l  l  ti  aict  0aB 
"  i=l  a=l  Z?«  j1eJa  j26JB  1  lct  aB 

*  n  X  gkk  gkk  (ffLL0ii 

n  k1,k2,k3Jk4  k1k2  k3k4  *1 *2  34 


*  *  ■  *  * 

+  &ili4  0i2X3 

(6.4) 

+  ai  i  TTi7ri+aii  11  i  77  i  +  a  i  i  71  i  77  i  +  °  i  i 

1  2  3  4  1  3  2  4  1  4  2  3  2  3  1  4 


+  CTi  i  7ri7ri+0ii  *i  *i  +  17  i  "i  TT  i  "i  ) 
2  4  1  3  3  4  1  2  1  2  3  4 


k  r  r 


3(t)  -  E^Ct)  —  III  l  l  t.  a. 

2»/n  1=1  a=l  6=1  jieJa  j2<-Jp  1  1 


rl  v  ,  *  *  *  * 

I  —  )  g,  ,  g,  .  ( a  a  +  a  a 

n  &k,k„  ^kok,  m.iru  num.  num«  m. 


k^kglkg,^  K1K2  k3K4  mlm2  ^4  mlm3  m2m4 


*  * 
+  a  a 

mlm4 


jft  jfc 

+  0  It  IT  +  0  IT  TT  +  0  IT  IT 

m^mg  nig  m^m3  m2  ^1^4  ^2  ^3 


+  0  tr  IT 
m2m3  m4 


j|c  jfe 

+  0  irir+a  tt  it  +  it  n  it  ) 

^2  m3  nig  rn ^  nig  ni^ 


u  y  g,  ,  (0  +  ir  it  )  -  a .  1 

j2  k1(k2  klk2  mlm2  ml  m2  J1  k  k4 


%k4  (om3m4l  V-4>  +  n 


where  ®  is  the  Kronecker  product,  and  as  in  (3.7) 


B 


1  ■  jlj  <i1aK(cU . clp> 


G  = 


~  =  ^Ul*  H2' * ' ' * Hn^ ’ 


(6.5) 


IT  = 


2i 


Cl  -  (I®  E)(G®  B1)l"1  K 

/n 


=  <VT2 . %>’ 


E~  =  (I®S)n  -  2-  fO®B.)(I®£)]-J-=  <ff*  ) 

Sn  A 


-1 


✓FT 


and 


ij  =  P  x  (kj-1)  +  J1 


m 


1  =  P  *  (k1-l)  + 


i2  =  p  x  (k3-l)  +  jj 


m2  =  p  x  (k2-l)  +  j 


io  =  P  x  (k0-l)  +  j. 


=  p  x  (ko-1)  +  j. 


ia  =  P  x  (k.-l)  +  j, 


m4  =  p  x  (k4-l)  +  j. 


Use  the  expansion  that 


d  -  2-  (I®E) (G®B..)1  x=  l  I  -1  (I®E)(G®B  )f 
/n  s=0  /n 


,-l 


2i 


i-4 


i  00  (2i)str  Gs •  tr(B.  E)s 

1 1  -  -  (G®B1E)|“  =  exp  4(1  - - - ) 

S=1  s^s 


(6.6) 


and  tr  W  (G®B1)  =  tr  B-^GU’  . 

We  obtain  the  characteristic  function  as 


f(t1,...,tk)  =  exp(-^  t’Qt) 
k 


(1  +  -  I  i  t.(h  +h2)  +  l  l 
11 


(i3t.  t  t  ) 
11  12  x3 
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where 


\ '  VN«,* 


(hg+h4+h5)l  +  0(n-i) 


(6.7) 


.  „2  ( i, )  ( i9)  ( *1 )  (2  i„) 

Q  =  (Qi  •  ),  Q,  t  =2  ££ L  tr  R  1  R  2  +  4  tr  R  1  2 

1112  1 i  2  n 


hl  =  l  l  l  a  0  J  (^-  (o  .  a,  +  a2  ) 
a=l  3^a  3  J2J2  J1^2 

J2rJ3 

+  2  \ia$x +  ^  °h,2  fjs 


h2  "  I  I  a  (~p-  a2  +  2  T(.2),  a  ) 

^  (X  3=1  -j  frj  -*-a”  n  ^1^2  JiJo 

>  *-*  i  “  cj  **■  “  i  « 


^2C  J8 


4  trG°  (i-i)  (i^)  (^3)  (i-i)  (ip)  (3,i3) 

h3  =  |  tr  R  1  R  2  R  3  +  4  tr  R  1  R  2  f  3 


(3.8) 


4  l  l 
a=l  8^a 


l  a.  0-1  (£*£.„'-*'+  fl'.T) 

d^Ja  Xla  aB  n  “0^2  ^2  Vl 


2  (i9)  (i„)  (i0) 


J2f  JB 


(ir2!5'Y+  4Y+  «-Y> 

^  J  ^  J  2  Jg^l 


2  y  y  (trG2  „(i2)  (12) 

A  L  aiaB  ^  n  ^11  *  ~  ^  1  i  ' 

a, 6  *af5  n  J1J1  J1J1 


^2eJ8 


trG2  ^3^  ( ^  3) 

(~n~  ~  j  ?  +  2  ni  ?  > 

n  j2j2  j2j2 


T  ,  R(1)  =  C(i)E,  C(i)  =  diag(ciJL>  .  .  .  ,cjp) 

(j>i)=  c(i)  OgV  f  H(i)  =  ER(i)f  0(i)  =  t(2)  R(i) 


(6.9) 


where  A.,  denotes  the  (i,j)th  elements  of  matrix  A  =  (A...) 

••■J 
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7.  Asymptotic  Distributions  of  Functions  of  the 
Roots  of  the  Complex  Wishart  Matrix 

Let  Z  =  +  i  Zg  be  a  p*n  matrix  and  let  the  rows 

cif  (Zj:  Zg)  be  distributed  independently  as  multivariate 
normal  with  covariance  matrix 


and  let  the  mean  vector  of  j-th  row  of  (ZJ :  Z' )  be 

I  I 

1  /  1  ^  (O')  -  - 

u.  =  (y.  *  ,  yv.  ).  Also  let  S  =  ZZ'  where  Z  denotes  the 
~J  ~J  ~3 

complex  conjugate  of  Z.  Then,  the  distribution  of  S  is 
known  to  be  central  or  noncentral  complex  Wishart  matrix 

with  n  degrees  of  freedom  according  as  M  =  0  or  M  =  0.  where. 

M  =  UW  ,  E(Z)  =  LI.  The  expected  value  of  S  is  given  by 

E(S)  =  2n( E^-iEg)  +  M 

The  matrix  S  is  Hermit ian  and  the  eigenvalues  of  S/n  are 
denoted  by 

^"3  —  ‘ 

In  the  sequel,  we  assume  that  E(S)  =  n  diag.  (A-,...,A  )  and 

x 

A. _>. .  . _>A  .  In  addition,  we  assume  that  A^s  have  multiplicity 
*  P 

as  in  (3.1).  Now,  let 

Lj  =  ^  tTj( ^i» • • • » £p)  “  Tj(^i . Ap)}  (7. 

for  j  =  l,2,...,k  and  the  function  TjC&j , . . . )  satisfy  the 
assumptions  (3.2)  and  (3.L).  Then,  following  the  same  lines 


as  in  Section  3  for  the  real  case,  we  obtain  the  following 
asymptotic  expression  for  the  joint  density  of  ,  L2>...,L.  : 


1  k 


^(L,.  •  •  •  ,Lfc)  =  N(L,Q)  [1+  -  l  Hi(L)(d1+d2) 

/h  i=l  x 


+  0( n_1 ) 


(7.2) 


where 


(i,)~(i9)  -  <i_ )-(ip) 

<*•  n  t  A  n  •*-  T5  " 


Q  ~  (Ql1i2),Qiii2  =  4  tr  Ri  1  Rj_  2  +  4  tr  R2  1  R 


(i1)-(i2) 


+  8  tr  R1  x  $ 

*v 

ana  Q  is  assumed  to  be  nonsingular.  Also, 


r  r 


d,  =  4  l  l  l  l  SS1}  ta'1’  i,  .  v  . 

1  i  a  *  -i  a  s  t  i®  JoJo  J-iJi  JoJp  3  9.3  9.  3 1  3 1  ' 


a=l  flfa  j1«Ja  J2pjc 


'lJl  J2J2  J1J1  J2J2  J2')2  Jl'il 


(7.3) 


^  r  r 

dp  =  2  I  l 


ll  l  l  a.  (ffi1!  '  +  2^i>  v  -  a\2).  ) 

ipIN  VJo  VJB  10B  ^  12  12 


(7.4) 


o  (I,).  (i„)„  (io)  .  ( i-i  )  (i2).(io) 

g1  =  f  tr  ft1  3  R3  2  Rx  J  +  8  tr  R1  -1  R^  *  d 

(I,),  do),  (ip)  .  ( i-i  ) .  (i2)-(i3) 

+  8  tr  R2  1  R2  2  R1  -  8  tr  R2  1  R2  *'  ' 


(7.5) 


.  r  r  ~(i2)  -(ip)  ~^2^  ~^i2\ 

go  =  is  y  y  y  y  a.  c(s,  t  +  t,  f  +  t .  f  +  a  .  ) 

g2  a=l  Bfa  j^Jg  ila  dlj2  *1J2  j2>H  -1J2 
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-^3)  -(*3)  -(^3)  -do) 

*  (=4  4  +  T.  .  +  T.  3.  +  G,  V 

^1^2  ^1^2  3  2^1  J 1 .1 2 

~(i2)  v(io)  .(iq) 

+  (Uj  j  ~  U1  1  )(Ui  i  "  Ui  i  >'l 
J1J2  J2J1  ^1^2  ^2^1 


2J1  J1J2 


(7.6) 


S3  ■  8  f  I  X  la,,  (f 

a=l  B-J  J2CJe 


-~(12>.  „r<V.  =«2> 


+  ^iV  °3131> 


(7.7) 


where  C 


(i) 


~d~)  -do)  .do) 

x-(-1  i  +  2T  3  +  G  ^  ) 

J2J2  J2J2  ^2^2 

,(2) 


=  dlaB  (CH . V’  ri-‘;ijia>.  Z2=(oi"i2) 


n 

I 

j=l 


M=  I  C^/^’+H/V2)')/-=^ij2) 


5  (i)  _  r( i). 


1  ~  ^  *  ^2 


(i)  =  C(l)r2>  tj)(i)  =  c(1) 


M 


f(i)  =  g(i)  =  ^  c(i)  Ef 


T(i)  =  m  c(i)t:  ,  u(i)  =  m  c(1)  j:, 

J-  < 


Krishnaiah  and  Lee  (1977)  derived  the  asymptotic  joint 
distributions  of  the  linear  combinations  and  ratios  of  the 
linear  combinations  of  the  eigenvalues  of  the  central  complex 
Wishart  matrix  when  the  roots  are  simple.  These  results  are 
special  cases  of  the  results  in  this  section. 
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